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Non-dissipative dynamics of interacting electrons in two tunnel-coupled quantum dots is studied the- 
Qs, ' oretically within the framework of the Hubbard model. Various values of intra-dot Coulomb repulsion 

0^ , energy U and inter-dot tunneling energy V are considered, which correspond to various size of the dots 

and to various distance between them. In the ground state, the average value of the spin projection (mag- 
netic moment) at each dot is zero. The input signal (the local external magnetic field H) applied to one of 
the dots at a time t — causes the electronic subsystem to evolve in such a way that magnetic moments 
of quantum dots become oriented in the opposite directions at any time t > 0. For any set of U and V, 
■ there exist optimal values of H and t which maximize the absolute values of magnetic moments at both 

dots, and magnetic moments become almost saturated. Thus, the antiferromagnetic-like spin ordering 
can be realized at the stage of coherent temporal evolution, well before the relaxation to a new ground 
state at the sacrifice of inelastic processes. This effect ("dynamical antiferromagnetism" ) may be used for 
implementation of a logic function NOT in an extremely short time. A possibility to use the arrays of 
C/2 \ quantum dots as high-speed single-electron devices of new generation is discussed. 
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1. Introduction 



Feynman's original idea [GJ to use the states of a quantum mechanical system for information coding and processing 
has stimulated the work on a quantum theory of computation 0. This theory takes advantage of fundamental 
properties of quantum systems, such as superposition, interference, entanglement, nonlocality, and uncertainty. The 
i i ' quantum bit ("qubit") can exist in arbitrary superpositions of classical bits, and 1. In this way it appears to be 
possible to develop very efficient computational algorithms which are exponentially faster than conventional ones 
H . The first experimental demonstrations of quantum gates (elementary units capable to perform simplest logical 
, operations) have been reported recently They are based on encoding the qubits in photon states. However, a 

lot of skepticism has been given to a possibility of practical utilization of such gates in a quantum computer [^,[7| . 
t-H , Another way is to use quantum mechanical systems to employ classical Boolean logic. For example, Lent et al. Q 
' have proposed to represent the bits and 1 by two opposite charge polarizations of a cell containing five quantum dots 
occupied by two electrons. Proper lateral arrangements of such cells with respect to each other allow the realization of 
various logical functions (NOT, AND, OR, etc). Nomoto et al. M have studied a logic device AND/OR composed of a 
pair of tunnel-coupled quantum dots. They viewed unoccupation/occupation of a quantum dot by a single electron as 
a bit 0/1. Bandyopadhyay et al. ||lo|| have discussed a concept of spin gates. In that approach the bits of information 
are carried by spins of individual electrons, e.g., the logical 1 (0) corresponds to the spin "up" ("down") direction 
1 ■ of electron spin at a given quantum dot. In arrays of quantum dots, the quantum tunneling of electrons between 
adjacent dots and/or Coulomb interaction of electrons with each other play a role of "wiring", resulting in the signal 
propagation from dot to dot ('quantum coupled architecture'). 
q \ Of course, the classical computational algorithms are not as efficient as the quantum ones since they do not benefit 
from purely quantum effects (superposition, interference, etc.). Nevertheless, one would expect that the reduction in 
. £h ! size of computer components will result in increase of the computational speed (up to inverse atomic times) and data 
storage density (up to several units per nm 2 for planar structures). Besides, the Boolean logic seems to be more closer 
to a practical implementation at a nanometer scale. Indeed, a number of techniques have been developed to produce 
arrays of quantum dots on a solid surface (see, e.g., 

In addition, the scanning tunneling microscopy provides 
a direct way to write-in and read-out the information since the tunnel current is sensitive to the magnetic moment of 
a single atom fill . Surely, the fabrication of classical logic devices based on quantum dots arrays will require more 
advanced technology. However it would be prudent to begin thinking now about the grounding in theory. 

Let us discuss more closely the spin gates proposed by Bandyopadhyay et al. and later investigated theoretically 
by Molotkov and Nazin [Q and by Krasheninnikov and Openov [jl5| . Such gates are assumed to consist of a number 
of tunnel-coupled quantum dots fabricated on a solid surface. The quantum dots play a role of potential wells for 
electrons. Each dot is supposed to have a single size-quantized energy level and there is, on the average, one electron 
per dot (the total number of electrons in the array of quantum dots can be controlled by adjusting the substrate 
voltage). In a spin gate, the stored information is determined by the spin configuration of the ground state of the 
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gate, i.e., by ground-state averages (S Z i) of spin projection operators S Z i at some dots (i is the number of a dot in 
the array). 

Each spin gate has input and output dots. The input dots serve for writing the information to the gate through 
appropriate orientation of their electron spins (e.g., by the action of the local magnetic field generated by a magnetic 
STM tip). The external influence changes the Hamiltonian of the gate. The new ground state (and its spin configu- 
ration) is different from the initial one since there are electron interactions within (Coulomb repulsion) and between 
(tunneling and Coulomb repulsion) dots in the gate. The values of (S Z i) at the output dots represent the result of 
calculations. The number of input and output dots in the gate depends on a particular logical function realized by 
the gate. For example, the AND gate has two inputs and one output, etc. 

The "ground state computing" is the key idea for implementation of classical Boolean logic on a nanometer scale 
(we note that this idea was used not only for spin gates, but also for other kinds of quantum-dots-based gates, see, e.g., 
IIP])- However, the concept of ground state computing has some essential drawbacks. The first one stems from the 
fact that the logical variables 1 and are associated not with certain pure quantum states, but with some quantum- 
mechanical ground-state averages. For a particular case of spin gates, the logical 1 (0) at the dot i corresponds to 
a positive (negative) value of (S X i) or, more precisely, to (S Z i) > St/2 ({S Z i) < —St/2), where St is the "threshold 
value", < St < 1 ]lj,|l5|. We note, however, that in a complex quantum-mechanical system, the electron spin 
at any quantum dot is never directed strictly "up" or "down", so the value of St cannot be taken equal to unity. 
Theoretical analysis of several logical spin gates (NOT, AND, OR, NAND, NOR, XOR, NXOR, and half-adder) within 
the framework of spin- 1/2 Heisenberg model has shown that for most of those gates the truth tables (which establish 
a one-to-one correspondence between the values of input and output bits) can be realized at sufficiently low values of 
St only (St = 0.1, or even St = 0.05 for half-adder) |Q. But the value of St should not be too low. Indeed, a low 
value of St does not imply that average projection of electron spin (S Z i) at the dot i has low absolute value (so that 
it is hardly seen by a reading device). In fact, since an electron spin is a quantum- mechanical entity, we can detect 
only its "up" or "down" directions at any measurement. A low absolute value of (S Z i) just means that a probability 
j>l to find an electron in the state with spin " up" doesn't differ much from a probability p± to find an electron in the 
state with spin "down". But if so, we will not be able to draw a definite conclusion about the sign of (S Z i) upon 
one particular measurement. In other words, there will be a high error probability P err , i.e., by a definition, a high 
probability to measure S Z i — —1/2 while (S Z i) > 0, or to measure S Z i = 1/2 while (S Z i) < (P err ~ 0.5 if St << 1). 
This is detrimental for the computational efficiency. 

For some simplest spin gates the issue of high error probability, related to low values of can be resolved 

by applying a suitably directed uniform external magnetic field to the whole gate C^). In particular, a fundamental 
possibility to reach ground-state averages | (S z i) | > 0.475 at input and output dots of NOT- AND and NOT-OR three- 
dot spin gates has been demonstrated theoretically in |Q. Such a large value of |(<S a <)| corresponds to P err < 0.05. 
However, apart from the problem of low | (S Z i) | , the concept of ground state computing has one more serious limitation. 

The point is that under the influence of external source on input dots of a spin gate, the relaxation of electronic 
subsystem to a new ground state can take place only through energy dissipation. The rate at which dissipation 
processes occur determines the speed of operation. However, whereas in two-dimensional nanostructures (quantum 
wells) the prevailing relaxation process is the emission of a longitudinal optical phonon with high relaxation rates of the 
order of 10 12 s" 1 , this process is forbidden in quantum dots because of the discrete nature of energy spectrum. Nomoto 
et al. H have shown that in semiconductor-based quantum dots, the emission of longitudinal acoustic phonon is the 
main contribution to the energy dissipation process. The relaxation rate is (10 6 -j- 10 9 ) s _1 , depending on geometrical 
parameters of an array of quantum dots P). The operation speed less than 10 9 s _1 is obviously too low for quantum- 
dots-based spin gates to be seriously considered as logic devices of new generation. Though some other scattering 
mechanisms (such as interface-phonon scattering) may also contribute to the energy dissipation in arrays of quantum 
dots H, it is unlikely that they will cause the relaxation rate to increase significantly. Hence, it is reasonable to search 
for some other principles of spin gates operation, beyond the concept of ground state computing. 

Recently Bandyopadhyay and Roychowdhury |l^| have studied theoretically a possibility of non-dissipative operation 
of the simplest two-dot spin gate NOT (inverter). They showed that under the influence of local external magnetic 
field, the logical function can be realized at the stage of coherent temporal evolution of electronic subsystem in 
time less than 10~ n s, i.e., well before the relaxation to a new ground state takes place at the sacrifice of inelastic 
interactions. Moreover, it was shown [ fl6| that there exists an optimal input signal energy to achieve a complete 
switching of the inverter, i.e., to realize spin configurations with saturated spin projections ((S Z i) = 1/2 or -1/2) at 
input and output dots. The results obtained in |l(| are very encouraging. However, the authors of Ref. [jl6| used the 
spin- 1/2 antiferromagnetic Heisenberg model to describe the correlated electrons in quantum dots constituting the 
inverter. As is known jL7|, the Heisenberg model is just a limiting case of the half-filled Hubbard model at U >> V, 
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where U is the intra-dot Coulomb repulsion energy and V is the inter-dot tunneling energy. Since it is not clear a 
priori to what extent the Heisenberg model can describe the actual experimental situation, it is instructive to make 
use of more realistic Hubbard model, without any constraints on the ratio U/V. 

In this paper we study the non-dissipative dynamics of the inverter within the framework of the Hubbard model 
(preliminary results have been presented in Jl8||). We explore a broad range of U and V values which physically 
correspond to various size of quantum dots and to various distance between them (the Heisenberg model as a limiting 
case at U >> V is also considered since it allows for a simple analytical solution). We show that complete switching 
of the inverter can be realized at unrealistic values of U /V — and U/V — oo only. Nevertheless, it turns out that 
at any value of U/V one can reach averages KS^)! > 0.45 at input and output dots in a very short operation time of 
the order of 1CP 11 s. Such large values of correspond to the error probability P err < 0.1 which is sufficiently 

low for implementation of logic function. 

The paper is organized as follows. The Hubbard and Heisenberg models for the two-dot spin gate NOT (inverter) 
are described in Sec. 2. The ground state computing is briefly outlined in Sec. 3 for completeness sake. The temporal 
evolution of electronic subsystem, once the local external magnetic field is applied to the input dot, is studied in 
Sec. 4, with special emphasis on time-dependent quantum-mechanical spin averages at input and output dots. We also 
analyze the dynamics of the inverter upon removing the input signal. The results obtained are discussed in Sec. 5 from 
the viewpoint of prospects for non-dissipative computing in spin gates based on quantum dots. In Sec. 6, concluding 
remarks are given. 



2. Theoretical models 

The inverter consists of two closely spaced quantum dots (A and B) occupied by two electrons M,0,E6|. One of 
two dots (say, the dot A) serves for writing the input signal to the gate by the action of the local external magnetic 
field Ha- The second dot (B) is the output. 



A. Hubbard model 

Under the assumption that each dot has only one size-quantized energy level, electron-dot interactions can be 
parametrized by one-electron Hartree-Fock on-site energies e^, es, and the energy V of electron tunneling between 
the dots . Electron-electron interactions are specified by the energies U aa , Ub b and U ab of Coulomb repulsion between 
two electrons residing at the same dot or at different dots. We ignore the exchange Coulomb interaction between the 
dots G3 since, first, its characteristic energy is generally much smaller than Uij (i,j — A, B) [|l9| and, second, the 
antiferromagnetic exchange naturally arises in the limiting case V < < Uij which corresponds to the reduction of the 
Hubbard model to the Heisenberg one [0 , see below. The presence of the (parallel to z-axis) local external magnetic 
field Ha at the dot A leads to the change in the total energy by [iHa , where \i is z-component of the magnetic moment 
at the dot A. In the site representation, the resulting Hubbard-like Hamiltonian has the form 

H = ( A flAa + £B ~ V J2( & Aa^Ba + & ta^Aa) + 

a cr a 

UAAnA]fl A l + U BB n B"\f^Bl + U AB ^ «Aa«B CT ' - 9^bH A S z a, (1) 

where a ia (af^) are the operators for annihilation (creation) of electron with the spin projection a = f or J, at the dot 
i = A or B, hi a = a~l a a i(7 are the electron number operators, S Z i — (h^ — n^)/2 is the operator of spin projection at 
the dot i, g is the Lande factor, fis is the Bohr magneton. The external magnetic field is directed along the magnetic 
moment of electron having a = f. 

Now let us slightly simplify the Hamiltonian (]!]), keeping all the key terms in it. First, we assume that the dots 
A and B are identical, so that ea — es = £o an d U aa — Ubb — Uq. Second, we suppose that the total number of 
electrons in the inverter N = n^t + f l A[ + + ^sj. equals to 2 (one electron per dot, i.e., the half-filled band) and 
remains unchanged during the operation at least till the reading of the result from the output dot. Then, by 

making use of the identity (see, e.g., |J) 

2 fiAafiBa' + ^n A ]n A i + 2h B ]h B i = N(N - 1), 

<T,a f 

we obtain 
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H = -V^2(a^ a a B(J + a+ a a Aa ) + (U - U A B){n A] n Ai + h B1 h Bl ) ~ 

a 

} d nAaSign(a) + e N + U AB , (2) 

a 

where sign(a) = +1 and -1 for a = j and J, respectively. Finally, omitting two last non-operator terms from the 
Hamiltonian (|2|), we arrive at the expression 

H = -V^2(a\ a a Ba + a+ a a Aa ) + U{n A] n Ai + h B] h Bi ) - H^n A aSign{a), (3) 

cr cr 

where we have introduced the renormalized energy of intra-dot Coulomb repulsion, U — Uq — Uab, and designated 
the input signal energy H = qhbHa/^- We note that generally U > since Uq > Uab @- 

The complete orthonormal set of inverter eigenstates is formed by the two-electron basis states 

|1} = | U>, |2) = | U>, |3> = | U,0>, |4> = |0,U>, |5) = | T,T>, |6) = | |,4>, (4) 

where, e.g., the notation | j, |) denotes the state with up-spin electron at the dot A and down-spin electron at the dot 
B, the notation | TJ.,0) denotes the state with two (up-spin and down-spin) electrons at the dot A and no electrons 
at the dot B, etc. We stress that the spin projection on either of the two quantum dots is zero in the basis states 
1 3) and |4). In this respect the NOT gate based on two quantum dots differs in general from the "spatially extended 
NOT gate" based on two spins |20| . In fact, the latter corresponds to the limiting case U » V of the Hamiltonian 
(||) and may be described within the framework of the Heisenberg model, see below. 

In the virgin state (i.e., at H = 0), the time-dependent electron wave function ^fo(t) is the solution of the Schrodinger 
equation 

ihd^ {t)/dt = H^ {t) (5) 
with the Hamiltonian (0) taken at H — and can be represented as 



*„(*) = E4 0) *i 0) expH4 0) <A), (6) 
fe=i 



where the coefficients Ar^ do not depend on time, and 'J'^ '' and (k = 1 -j- 6) are eigenvectors and eigenenergies 
of the stationary Schrodinger equation 



^ 0) = 4°M 0) (7) 



They are as follows: 



*i 0) = ^(|l)-|2)), ^ = -L(|3) - |4)), 



= - Ji + = id + 12) + ^+m>-u VOTM/ 



2 V 1 1 > + |2> 4F |3> 4F 1 



*<°> = |5), ^ 0) = |6), (8) 



and 



£f = 0, = C/, tff = -^ 7? +f 7 "- t/ , = ^^+^ = 0, = 0. (9) 

Generally, the coefficients A^ in (^|) can take any values which satisfy the normalization condition 

Ei4 0) i 2 = i- (io) 

k=l 
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In other words, the electronic subsystem can exist in a rather complex superposition of quantum eigenstates "J^ ' . 

In the presence of the input signal at the dot A, the time-dependent wave function ^(t) satisfies the Schrodinger 
equation (||) with the Hamiltonian (||), and can be expanded as 

6 

= A k * k exp(-iE k t/h), (11) 
fe=l 

where ^ k and E k (k = 1 — 6) are eigenvectors and eigenenergies of the stationary Schrodinger equation 

HV k = E k V k , (12) 
and the time-independent coefficients A k satisfy the normalization condition 

6 

£iV = i- ( 13 ) 

k=l 

At arbitrary values of U, V, and H the eigenvalue equation ([l2|) reduces to the algebraic equation of the third 
power in E k . The resulting expressions are too cumbersome for analysis, so it is more convenient to solve the equation 
( |l2| ) numerically. However, it is useful for the following discussion to consider a special case U = which physically 
corresponds to U << V (closely-spaced large-sized dots ||). At U = and at arbitrary values of V and H we have 
for eigenvectors and eigenenergies of the equation ( |l2] ) rather simple analytical expressions: 







|4>) 




_ 2V 


\Aff 2 + 


AV 2 - 


"*! 


V# 2 + 


AV 2 - 


fif 


\Aff 2 + 


AV 2 - 





V# 2 + AV 2 + H . . VH 2 +AV 2 -H. n , VH 2 + AV 2 -H .„. V# 2 + 4F 2 

*3 = — -===^- |1> + = : _ Q ; - 12) + — 13) 




2VH 2 + AV 2 \ VH 2 +AV 2 + H 2V' 1 2V 

_ y<H 2 +AV 2 -H ( y<H 2 +AV 2 + H VJT 2 + AV 2 + H y/H 2 + AV 2 + H \ 
4 ~ 2VH 2 + AV 2 \ ' + VH 2 +AV 2 -H l ' 2V l) 2V ''J' 

*5 = |5), * 6 = |6>, (14) 



and 



E x = 0, E 2 = 0, E 3 = -y/H 2 + AV 2 , E± = y/H 2 + AV 2 , E 5 = -H, E 6 = H. (15) 



One can see that expressions (|lj) and (|l5|) taken at H = coincide with expressions (H) and (Q) taken at U = 
respectively . 



B. Heisenberg model 

In the limiting case U » V (widely-spaced small-sized dots H), the half-filled Hubbard model reduces to the 
spin- 1/2 antiferromagnetic Heisenberg model (l7| . The corresponding Hamiltonian for the inverter has the form 

H= o-gA^aB - Hy^h Aa sign{a), (16) 

where a a i (ct = x, y, z) are the Pauli matrices describing the electron at the dot i = A or B and J = V 2 /U is 
the energy of exchange interaction between two electrons residing at dots A and B. We have omitted the term 
— J(fiAi + n,Ai)(nB\ + ft-Bi) 0] from ( |l6| ) since in the Heisenberg model the complete orthonormal set of inverter 
eigenstates doesn't include the two-electron states |3) = | fl, 0) and |4) = |0, t j), see (||), so that each dot is occupied 
by one electron, and hence (riytf + fiAi){nB^ + Ubi) = 1 is just a c-number. 

The Heisenberg model for the inverter allows for analytical solution at any values of J and H @,^6| (in what follows 
we shall make use of this fact to analize the case U >> V). The eigenvectors and eigenenergies of the stationary 
Schrodinger equation with the Hamiltonian (]l^) are as follows: 
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* 3 = |5), * 4 = |6), (17) 



and 



Ex = -J- V^ + 4J2, E 2 = -J+ ^H 2 +4 J 2 , E 3 = J-H, E 4 = J + H. (18) 



The time-dependent wave function \I/(t) has the form (jll]) with 'J'fc and i7fe given by ( |l7| ) and ( |l8| ) respectively 
except that the number of inverter eigenstates in the Heisenberg model is less by two than in the Hubbard one, so 
that the summation over k in ( [ll] ) is from k = 1 up to k = 4. 

3. Ground state computing 

The truth table of the logical gate NOT is very simple, the output bit B being just the inverse of the input bit A, 
i.e., B~A(B = l\iA = Q and B = if A = 1). Let us see how this truth table is realized if the input and output 
bits are encoded in ground state averages (S z a) and (S z b) respectively. We shall ascribe the bit 1 (0) at the dot i to 
the positive (negative) value of (S z i), with the constraint > St/2, where < St < 1 

As discussed in the Introduction, one wishes to have St as large as possible since the greater is St, the lower is the 
error probability P err , i.e., the probability to read the "wrong" signal S Z i at the dot i (e.g., to measure S z i = — 1/2 
while (S z i) > 0). The physical reason for a possibility of such an error in the device operation is the entangled structure 
of the ground state wave function (the latter cannot be factored as | |)s or | [)a\ T)b)- It is straightforward to 
find the upper limit of P err as a function of St- For the sake of definiteness, let us consider the case (S z a) > St/2, 
(S z b) < —St/2 which corresponds to H > 0. Then P err = 1 — p-f j, where p^ is the quantum-mechanical probability to 
find up-spin electron at the dot A and down-spin electron at the dot B. Since, by a definition, (S z a) — (j>n —Pii)/2, 
one has generally p*n — > S t , and hence P err < 1 — S t (in the special case U » V, the ground state eigenvector 
is composed of states | T,J.) and | J.,f) only, see (f|), fll~7j), (|l8|), and one has p^ + Pii = 1, Pu > (1 + S t )/2, and 
Perr < (1 — St)/2). Note that P err = in the ideal case St — 1 only, i.e., in the case of saturated spin projections 
\{S zi } \ = 1/2 at both dots i = A and B. 

It follows from (|J) and ([)]) that at zero input signal H — the ground state of the inverter at any value of the 
ratio U/V is the entangled state ^E'g '' with zero magnetic moments at both dots A and B since (^^ISzaI^^) = 
(tyf^lSzBl^f^) =0. At H 7^ 0, the ground state averages of spin projections at the input and output dots become 
nonzero. They are equal in magnitude and opposite in sign, (S z a) = —(S z b), just as required by the truth table of 
the inverter. Making use of d|) and (||), one has 

(§ X J = —-JL==, U«V, (19) 



2VH 2 + W 2 ' 



while 



TT 

(Sza) = , U»V (20) 

V 1 2VH 2 + 4 J 2 V ' 

as follows from ( p~7|) and (|l^). At arbitrary value of the ratio U /V, the dependence of (S z a) on H may be approximated 
as (S zA ) = H/2^/H 2 + m 2 where W = Vf(U/V), and the function f(x) is shown in Fig. 1 (/(0) = 1; f(x) = l/x 
at x » 1). Thus, (S z a) increases monotonically with H from (S z a) = at H = and asimptotically approaches its 
maximum value (S z a) = 1/2 at H >> W. The error probability is close to 0.75 (for U « V) or 0.5 (for U >> V) at 
H « W and decreases with H as P err ~ aW 2 /H 2 at H » W (a = 2 and 1 for U « V and U » V respectively). 

In order to minimize P err one should have the input signal energy as high as possible. However this cannot be 
easily realized experimentally since realistic values of Ha < 1 T correspond to quite low values of H < 0.1 meV, 
while H should be much greater than the characteristic energy W. On the other hand, the energy W should exceed 
considerably the thermal energy ksT for stable operation so the operating temperature should be restricted to 
at least 100 mK. 
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Aside from the problem of P err reduction, the speed of device operation is limited by the finite rate of relaxation of 
electronic subsystem to a new ground state upon applying the input signal (see the Introduction) . The relaxation rate 
is determined by dissipation processes, so the dissipation is a necessary attribute of the ground state computing in 
the array of quantum dots. In what follows we study a possibility to perform a logical operation NOT in the absence 
of dissipation, at the stage of non-dissipative temporal evolution of electronic subsystem. 

4. Non-dissipative computing 

Our purpose is to calculate the quantum-mechanical averages 

S zA {t) = S zB {t) = (*(t)|& A |*(t)) (21) 

at t > 0, i.e., upon applying the input signal to the dot A. For definiteness, we shall consider the case H A > (i.e., 
H > 0). We shall also assume that at t < (in the absence of external signal) the system is in its ground state, so 
that the time-dependent wave function at t < is <&o(*) = \tg exp(— t/K), see (^), (^|), and hence the coefficients 
A { ° ] in the expansion of ¥ (i) (§) have the values A { ° } = A^ = A ( ° ] = A { ° ] = A { ° ] = 0, A { ° ] = 1. 

The coefficients Ak in the expansion of the time-dependent wave function \&(t) ([ll]) at t > should be found from 
the initial condition = 0) = ^o(t = 0). It is convenient to write the eigenvectors of the stationary Schrodinger 
equation (|l2|) as 

6 

*fc = H B fc»- ( 22 ) 

n=l 

Then 

6 

71=1 

where 

6 

fn(t) = A k B kn exp(-iE k t/ti). (24) 
fc=i 

The probability to find the system in the basis state |n) at time t is 

Pn(t) = \fn(t)\ 2 . (25) 

Since, as follows from (|), (l\S zA \l) = 1/2, (2\S zA \2) = -1/2, (3|& A |3) = (4|5,a|4) = 0, (5|&a|5) = 1/2, (6|S Zj4 |6) = 
-1/2, one has from @, (§|), (§§): 

o zA {t) = -S zB [t) = . (26) 

First we consider the limiting case U << V |lj] (the weak coupling limit of the Hubbard model). Setting [7 = 0, 
we have rather simple equations for the probabilities p n (t) : 

■sin 2 (wt/2) 



Pl(t) 


40 


P2(t) 




P3(*) 


= P*(«) 







ij2 + 4^2 

1 / 16iT 2 V 2 . 4 , 



(27) 



where 



y/H 2 + 4V 2 /h. (28) 



7 



From (E6j) and (27) it is straightforward to find that 



S zA (t) = -S xB (t) = H l H + \ V 2 sin 2 M/2). (29) 



It follows from (29) that at any positive values of H and V the function S X A(t) is non-negative and oscillates in 
time with the frequency lo (given by (^)) and the amplitude 

SMto) = H l H + V w2 , (30) 

where to stands for the time of the first peak of S z A(t), 

7T Trh 

to = - = , (31) 
w VH 2 + W 2 

The curve S Z A(t) is shown in Fig. 2 for several different values of H/V. From Fig. 2 and fl2S|), ( |3l| ) one can see that 
the oscillation frequency lo monotonically increases (i.e., the time to decreases) with H, while the amplitude S Z A(to) 
first increases with H, reaches the maximum value 1/2 at H = 2V and then decreases, being in inverse proportion 
to H at H » V. Thus, a complete switching, 5 z a(*o) = 1/2 and 5 z s(to) = — 1/2, is achieved at H/V = 2 and 
t = nh/2V2V. 

It should be stressed that though the function S z A(t), being periodic in time, has its peak value at tk = 7r/w + 27rfc/a; 
(k is an integer), here and below we are interested in the lowest possible value of tk (i.e., to) since we want not only 
to reach the maximum permissible value of S z a, but to do it in as short as possible switching time. 

Now let us turn to the opposite limiting case, U >> V (the strong coupling limit of the Hubbard model). This case 
may be described within the framework of the Heisenberg model (Q) . It has been considered by Bandyopadhyay and 
Roychowdhury in flq| . For the sake of completeness, we briefly sketch the results obtained in |fi6[ . 

As we have noted above, in the Heisenberg model, the complete set of inverter eigenstates doesn't include the states 
|3) = | Tl,0) and |4) = |0, ||), see (f|), (|l7|). For the probabilities p n (t) to find the system in the basis states \n) = 
jl), |2), |5) and |6) one has 

Pi(t) = - 1H 5 r sin 2 ( ut 2) 

y W 2 V H 2 +4J 2 1/7 

P^) = -{l- w ^ J - 2 ^t/2) 

P 5 (t) = Pe(*) = , (32) 



where 



lo = 2VH 2 +4J 2 /h (33) 
(we recall that J — V 2 /U <<Vis the energy of antiferromagnetic exchange). From (pq) and (p2) one has 



S zA (t) = S zB (t) = H 1 H + J AJ2 sinVt/2). (34) 

One can see from ( ]3~i| ) that the function S zA {t) at [/ >> V has the same form as in the case U « V (^§). The 
difference between the two cases is that V in ( ps|) and ( p9|) i s replaced by J in ( |33"1) and (Q) respectively, and that 
there is an extra factor of 2 in the oscillation frequency in (|33j) as compared with (|28|) . The amplitude of the function 

S zA {t) is 

S Mto) = H 2+ J 4J 2 ' ( 35 ) 

where to is the time of the first peak of S z A(t), 

TT Tift 

to = - = , (36) 
lo 2VH 2 + 4 J 2 V ' 

It follows from p^ ) and ( ^6| ) that a complete switching, S z A(to) = 1/2 and S Z B{to) = — l/2, is achieved at ff/J = 2 
and t = 7r?i/4V2J- 
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Thus, having considered two limiting cases of weak (U << V) and strong (U >> V) electron coupling, we see 
that the complete switching of the inverter can take place in either limit of the Hubbard model providing a proper 
choice of the values of the input signal energy H and the switching time to- Hence, one may expect that the complete 
switching can occur also in the case of intermediate coupling, at an arbitrary ratio of U to V (i.e., at arbitrary values 
of the size of quantum dots and the distance between them). 

To check this hypothesis, we have calculated numerically the dependencies of S z a on t at different values of U /V 
and H/V . Several typical examples of S z a versus t curves are shown in Fig. 3. The switching time to was generally 
defined as a time of the first peak on the curve S z A(t), to being a function of U/V and H/V. As opposed to the 
cases of weak and strong coupling considered above, at arbitrary values of U/V and H/V the function S z A(t) is not 
periodic in time since it includes several harmonics with different frequencies and amplitudes. Hence, in principle, 
the value of S z A(t) greater than S z A(to) can be achieved at some longer time (see, e.g., Fig. 3b). But this case is of 
no interest for us since our purpose is not only to maximize S z a, but also to minimize the switching time. 

The curves of S Z A(to) versus H/V are shown in Fig. 4 for several values of U/V. One can see that an increase in 
U/V first results in a decreased height of the maximum S^ x on the S z A{to) versus H/V curve. For U/V > 2 the 
value of S^ A ax increases once more, but doesn't reach the saturation value 1/2 at finite U/V, though S^ A ax — > 1/2 if 
U/V — ► oo (this corresponds to the Heisenberg model and agrees with the results obtained in fl6[| ). Fig. 5 shows the 
dependence of S™£ x on U/V. Note that S™£ x > 0.45 at any U/V. 

We define H opt as the "optimal" value of input signal energy H, i.e., as the value of H at which the peak value 
S zA (t ) reaches its maximum S^ x at a given U/V (e.g., H opt = 2V at U/V = 0, see ©; H opt = 2J = 2V 2 /U 
at U » V, see ©; H opt = 0.83V at U/V = 2, see Fig. 4). The dependence of H opt /V on U/V is plotted in 
Fig. 6. One can see from this figure that H opt /V decreases monotonically with increasing U/V. It is remarkable 
that the equation H opt /V — 2V/U derived analytically in the Heisenberg model, i.e., for U >> V (p5|), is a good 
approximation of the numerically calculated curve even at U ~ V. It follows from Fig. 6 that the curves S z A(t) 
shown in Figs. 3a and 3b are for the case H > H opt while the curve S z A(t) shown in Fig. 3c is for the case H = H opt 
(hence, the peak value of S z A{to) in Fig. 3c is the maximum permissible value of S z A(t) at U/V = 1.5). 

An important characteristic of the logical gate is its switching time to which depends on U/V and H/V (see above). 
At a given U/V, we define the "optimal" switching time t opt as the time of the first peak on the curve S z A(t) at 
H = H opt . Hence, S z A{t op t) equals to S™ x , the maximum permissible value of S z a (to) at fixed U/V (see Fig. 4). 
The dependence of t op t on U /V is shown in Fig. 7. It is seen that t opt increases monotonically with increasing U /V 
from t op t = 7r?i/2v / 2V A at U/V = 0, sec (pH]), and approaches the curve t opt — nTiU / %\[2V 2 , derived analytically in the 
Heisenberg model, for U » V, see fl36|). 

It is instructive to examine the dynamics of the inverter upon removing the external signal from the input dot A at 
a time t — T without reading the spin polarization in the output dot B. At t > T the electronic subsystem continues 
to evolve in accordance with the Schrodinger equation (^) . For the sake of simplicity, we restrict our consideration to 
the case U » V (Heisenberg model) since this allows for a simple analytical solution. For the quantum- mechanically 
averaged spin projections at the dots A and B we have 

S zA (t) = S zB (t) = S sin ( AJ ^- T ) + arc tan ( —^L== tan(wT/2)^ , (37) 



VH 2 + 4 J 2 



where 



So = H 2 + 4J2 V / (^ 2 cos 2 (^T/2) + 4j2)sin 2 ( W r/2) (38) 

and uj is given by Eq. ([33|). We see that S z A(t) and S Z B(t) oscillate between — So and So having zero average values. 
The amplitude of oscillations So depends not only on the values of H and J but on the value of T as well. If H > 2 J 
then So has its maximum value 1/2 at uT = arccos(— AJ 2 / H 2 ) + 2irk where k is an integer. It is interesting that 
So = at loT = 2nk, i.e., S z A{t) = — S z s(t) = 0, just as in the ground state at H = 0. In this case the result of two 
subsequent perturbations at t = and t = T is such that the system at t > T behaves as if it was not perturbed at 
all. This effect can be explained in a way that the two perturbations interfere in a countcrphase thus compensating 
each other. 



5. Discussion 

In Sec. 3 we saw that in the presence of input signal at the dot A the ground state average of the spin projection 
(S z a) increases monotonically with increasing the input signal energy H from (S z a) = at H = and asymptotically 
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approaches the saturation value (S z a) — 1/2 at H >> W where W = Vf(U/V) is a characteristic energy which 
depends on U and V (W = V at U/V = and W = J — V 2 /U at U/V » 1, see Fig. 1). In the case of temporal 
evolution of electronic subsystem upon application of the input signal to the dot A one would expect the peak value 
S z A(to) of the time-dependent spin projection S z A(t) = (^(t^SzA^it)) be also governed by the same parameter, the 
ratio H/W. As follows from analytical expressions ( J30| ) and ( p5| ) obtained for S z A{to) in the limiting cases U/V « 1 
and U/V » 1 respectively, this is so indeed. However, unlike the ground state average (S z a), the dependence of 
S z A(to) on H/W is nonmonotonic. By contrast, S z a(£q) has a maximum S^ x at a finite value of H/W while tending 
to zero at H/W » 1, see (|30|), (|35|), Fig. 4, and Fig. 6. Hence there is a fundamental difference between the physical 
processes responsible for an increase of the spin projection in the static and dynamic cases. 

We stress that perfect antiferromagnetic spin ordering with maximum absolute values of spin projections on both 
dots, S z a = —S z b = 1/2, can be realized at the stage of coherent temporal evolution not only in the limit U/V — > oo 
(this would seem natural by analogy with the static case when large values of U/V are favorable for antiferromagnetism 
|p7| ) but also in the limit U/V — > 0. This effect, which we have called "dynamical antiferromagnetism", arises from 
the peculiar nature of electronic transitions under the influence of external perturbation. 

So, in the case of temporal non-dissipative evolution of the inverter, the maximum permissible (for given U and V) 
value of S z a is achieved at a finite value of the input signal energy H, i.e., S z a = S™j/ X at H = H opt and t = t opt , see 
Figs. 2 - 4. It is important that S™£ x > 0.45 at any ratio U/V while S™£ x -> 1/2 at U/V -> or V/U -> 0, see Fig. 
5. Since S Z B(t) = —S z A(t) at any time, then if the spin polarization (logical bit) at the output dot B is measured 
at a time t = t optl the inversed input signal, i.e., S z b = —1/2 for H > or S z b = 1/2 for H < 0, is detected with 
high probability. In other words, the error probability P err , i.e., the probability to measure the wrong output signal 
(S z b = 1/2 for H > or S z b = —1/2 for if < 0) is very low. Indeed, the upper limit of P err may be estimated as 
P err (t) < 1 — 2S z A(t), see Sec. 3. At t — t opt one has Perr{t op t) < 1 — 2S™j/ x , i.e., P e rr(t pt) « 1 at U « V or 
U >> V, the maximum value of P e rr{t pt) being about 0.1 at U /V = 2. 

We recall (see Sec. 3) that in order to reduce P err in ground state computing, one should increase the input signal 
energy H. In the case of interest, P err << 1, it can be shown from ([l9]) and ( p0| ) that H w W/^/P e rri where W is 
a characteristic energy that determines the scale of the energy levels scheme of the inverter, see Fig. 1. But if so, 
the value of H should be also much greater than the thermal energy ksT since the inequality ksT << W must be 
fulfilled in order the inverter was in its ground state. Thus, in fact, one should have 

k B T «W « H (39) 

in ground state computing. This condition is obviously hardly realizable since even the temperature as low as T = 1 
K corresponds to rather high magnetic field about 1 T, and hence at T = 1 K two strong inequalities in (|3^) can be 
satisfied only for unrealistic magnetic field about 100 T. Otherwise the error probability P err is very high. 

In the case of non-dissipative computing, there is no need for the relaxation to the new ground state since the perfect 
antiferromagnetic ordering of magnetic moments is realized at the stage of unitary evolution when an appropriate 
superposition of the inverter eigenstates takes place, see above. Hence, the input signal energy H should not be 
anomalously high. For example, in the limiting case U >> V, the optimal value of H is H op t = 2V 2 /U, so that, e.g., 
at V « 1 meV and U « 10 meV the external magnetic field acting on the input dot can be relatively low, about 1 T. 
Such a magnetic field is generated, e.g., by the magnetic moment of the order of hb in its vicinity. Hence, the input 
signal may be thought of as being supplied by the output quantum dot of another gate whose magnetic moment has 
been oriented in such way or another at the previous stage of calculations. 

We note that the limiting case U » V is much more suitable for the physical implementation than the case 
U << V. Indeed, at U << V we have H opt = 2V, hence for a realistic value of the input magnetic field about 1 
T, the Coulomb repulsion energy U should not exceed 0.01 meV in order to satisfy the inequality U « V . Since 
U ~ e 2 /ea, where a is the dot size and e is the dielectric constant J2l]], for U « 0.01 meV one has a w 10 4 - 10 5 nm. 
Such a relatively large structure obviously can not be considered as a quantum dot. 

Though the switching of the inverter to the new ground state upon the influence of the external signal is not 
necessary in non-dissipative computing, we stress that our results on unitary evolution of electronic subsystem of the 
inverter are obtained under the assumption that at t < (in the absence of the external signal) the inverter is in its 
ground state. Thus, one faces the problem of preparing such an initial state of the inverter, and consequently, the 
problem of temperature reduction, as in the case of ground state computing. In this respect, it would be interesting 
to search for such a type of input signal (e.g., for such a direction of the external magnetic field, see Mj) whose 
influence on the input dot will result in the desired temporal evolution of electronic subsystem (and consequently, in 
the realization of the truth table) irrespective of the initial state of the inverter. 

The other facing problem is the reading of the result of computation (output). While in ground state computing 
the output can be read at any moment after relaxation of the system to its new ground state, in non-dissipative 
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computing the reading time must be chosen with a great accuracy as a time when spin projections at the dots have 
their maximal values, in order to minimize the error probability. It is interesting to analyze the relationship between 
the "optimal" values of t opt and H opt . The dependence of the product t opt H opt on U/V is shown in Fig. 8. One can 



see that the value of t opt H opt decreases monotonically with U/V from Tvh/y/2 at U/V = 0, see (|30|), (31) down to 
irh/2\/2 at U/V — + oo, see (|3^). Thus, the relationship 

t opt « h/H opt (40) 

holds true for the whole range of U /V. In general, this is a consequence of the Heisenberg's uncertainty principle. 

We recall that our purpose is to minimize the switching time t opt and thereby to maximize the speed of computation. 
But according to j40|), the decrease in topt can be achieved only at the sacrifice of increase in H op t. Since the value 
of the input signal energy is finite, the relationship imposes a restriction on topt- For realistic value of external 
magnetic field ~1T one has t op t ~ 10~ n s. 



6. Conclusions 

The system consisting of two tunnel-coupled quantum dots, input and output ones, occupied by two electrons has 
been considered within the framework of the Hubbard model with arbitrary values of the parameters, the intra-dot 
Coulomb energy U and the inter-dot tunneling energy V. After applying the external signal (the local magnetic field) 
to the input dot, the magnetic moments of the dots become oriented in opposite directions, i.e., align antiferromag- 
netically. This allows one to use such a system as the quantum spin gate " inverter" realizing the truth table of logical 
NOT. 

In ground state computing, the inverter switches (relaxes) to the new ground state under the influence of the input 
signal. Since the relaxation process is due to inelastic interactions, the switching time is of the order of 10~ 9 -r 1CP 6 s 
which is rather long on the atomic scale. Moreover, the magnetic moments of the dots are far from being saturated, 
giving rise to the large probability of reading the wrong result of computation from the output dot. For these reasons 
the principle of ground state computing seems to be ineffective for the spin logical gates. 

We have suggested a new approach which overcomes the two above mentioned drawbacks of ground state computing. 
Our approach takes advantage of the coherent non-dissipative evolution of electronic subsystem after applying the 
signal to the input dot at t — 0. We have demonstrated that at any moment t > magnetic moments of the dots 
arrange in antiferromagnetic order, i.e., S z A(t) = —S Z B(t), thus realizing the truth table of the inverter. The value of 
SzA^t) (and hence S Z B(t) as well) is a nonmonotonic function of time, being a periodic function in the limiting cases 
U << V and U >> V. At given values of V, U, and the input signal energy H there always exists a moment to which 
corresponds to the first local maximum on the curve S z A(t). Besides, there exists an "optimal" value H opt at which 
the magnitude of S z A(fo) reaches its maximum permissible (for a given ratio of U/V) value S^ x , thus reducing the 
error probability P er r- Our calculations showed that S™ A ax > 0.45 at any U and V, being close to 1/2 at U << V 
and U » V . So, the value of P err << 1 is well below than that in ground state computing. The "optimal" values 
of H opt and t opt — to(H opt ) are related by Eq. (f40|). For reasonable values of U, V and realistic magnetic fields ~ 1 
T one has t opt ~ 10~ n s, several orders of magnitude less than in ground state computing. 

Note that the effect of "dynamical antiferromagnetism" appears to be usable for further development of more 
complicated spin gates consisting of a greater number of quantum dots. However, several problems should be resolved. 
First, our results are obtained for a particular case that the initial state of the inverter is its ground state. It would 
be desirable to search for such a type of the input signal (e.g., for such an orientation of the external magnetic field 
with respect to the quantization axis) that the complete switching of the inverter took place irrespective of its initial 
state (the latter may be a superposition of inverter eigenstates). Second, there is a problem of reading of the result 
of computation. If the direction of the magnetic moment on the output dot is measured at a time different from the 
"optimal" time t op t, then the error probability increases. We note that the quantum-mechanical average of the spin 
projection at the output dot of the inverter has the same sign at any time after the input signal had been applied to 
the input dot (this sign determines the value of the logical bit at the output). Hence, should the electron spin be a 
classical vector, then the result of calculation could be determined unambiguously. It would be very interesting to 
search for such a measurement procedure that enabled us to draw a conclusion about the average direction of the 
magnetic moment over some period of time. Solving of these problems appears to be the target of future research. 
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Figure captions 

Fig. l. The function f(x ) entering into the expression for the ground state average of the spin projection (S z a) = 
H/2 V / H 2 +W 2 f 2 (U/V) at the quantum dot A. 

Fig. 2. Temporal evolution of the quantum-mechanically averaged spin projection S z A{t) = (t)\ S z a \^{t)) at the 
quantum dot A upon applying the local external magnetic field Ha to the dot A. The intradot Coulomb repulsion 
energy U = at both dots A and B (the weak coupling limit of the Hubbard model). H = V (dashed line), H = 2V 
(solid line), H = 6V (dotted line). Here H = g^BHA/'i is the input signal energy, V is the energy of electron tunneling 
between the quantum dots, the time t is measured in units of fi/V . 

Fig.3. The same as in Fig.2, for a) U/V = 2, H/V = 2; b) U/V = 1, H/V = 4; c) U/V = 1.5, H/V = 1.155. 

Fig. 4. The peak value of S z A(to) versus H/V at different U/V. The switching time to is defined as a time of the 
first peak at the curve S z A(t) at given values of U/V and H/V. 

Fig. 5. The dependence of S^ x on U/V. Here S™j/ X is the maximum value of S z A(to) on the S z A(to) versus H/V 
curve at a given U/V. 

Fig. 6. The "optimal" value of normalized input signal energy H opt /V versus U/V. Here H opt is the value of H 
at which the peak value S'z^o) reaches its maximum S^ x at a given U/V, see Fig. 4. Solid line is the result of 
numerical calculations. Dashed line is the curve H opt /V = 2V/U obtained in the Heisenberg model, i.e., for U » V. 

Fig. 7. The "optimal" switching time t op t measured in units of Ti/V as a function of U /V . Here t op t is the time of the 
first peak on the curve S z A{t) at H = H opt , i.e., t opt (U/V) = to(U/V,H opt /V) corresponds to S™^, the maximum 
permissible value of S z A(to) at a given U /V. Solid line is the result of numerical calculations. Dashed line is the curve 
toptV/h = ttU/4:^/2V obtained in the Heisenberg model, i.e., for U >> V. 

Fig. 8. The product of "optimal" switching time t opt measured in units of h/V by the "optimal" value of normalized 
input signal energy H opt /V versus U/V. 
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